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Abstract 

O , 

' The ambiguity involved in the definition of effective-mass Hamiltonians for 

^SJ . nonrelativistic models is resolved using the Dirac equation. The multistep 

i approximation is extended for relativistic cases allowing the treatment of ar- 

bitrary potential and effective-mass profiles without ordering problems. On 
^ ' the other hand, if the Schrodinger equation is supposed to be used, our rel- 

ly-^ ' ativistic approach demonstrate that both results are coincidents if the Ben- 

. Daniel and Duke prescription for the kinetic-energy operator is implemented. 

Applications for semiconductor heterostructures are discussed. 
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The effective-mass theory has been successfully used in semiconductor heterostructures 
Pi I lU. An interesting aspect arises when we treat materials whose properties change from re- 

2 ■ gion to region. In particular, when the effective mass depends on position, the Schrodinger 
equation for an arbitrary potential profile is usually solved numerically by different methods. 
However, one of the problems of the effective-mass theory for semiconductors heterostruc- 
\ tures, is to decide how to write out the Hamiltonian operator. This problem arise from 
\ canonical quantization of the classical Hamiltonian. For position-dependent carrier effective 
mass, we have an ordering problem with the kinetic energy operator (KEO). Some authors 
proposed different forms for the kinetic energy operator, all having the generic form proposed 
by von Roos [H 
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- [nf" {x)pm^ pro' [x) + m? {x)pm^pm"{x)"^ , (1) 



where a + jS + j = —1, but the problem was not resolved because there is not a first principle 
to fix only one operator. 

This ambiguity indicates that the Schrodinger equation is not rigorously suitable in the 
effective-mass aproximation with position-dependent effective mass. It is reasonable to try 
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other equation that represents the same physics of the Schrodinger equation in the low 
energy hmit. 

The object of this Letter is to demonstrate that, in fact, the Dirac equation (at adequate 
hmits) can successfully be used to describe quantum-mechanical systems where position- 
dependent effective mass is present. We recall that, in Dirac equation, the kinetic energy 
operator and the mass term appear separately, so there is no ordering problems in this 
context. 

Obviously the considerations of this work only concern the mathematical issues related 
with the equations of motion. Indeed, a physically sensible application of the Dirac equation 
to semiconductor heterostructures would have to take in to account a relativistic extension 
of the Wannier-Slater theorem . 

In this work we use a numerical method (multistep potential approximation [Q) which 
has been applied to solve Schrodinger equation for an arbitrary potential profile. Here we 
extend this algorithm to the relativistic case, in such a way that the ambiguity problem is 
overcome. In particular, we consider the Dirac equation with a one dimensional arbitrary 
potential well and find the energy levels for a particle. Also we apply the method for a 
particular type of heterostructure and compare the results to those obtained in the context 
of the Schrodinger equation with the form (|ip for the KEO and several choices for the 
parameters and 7. For a KEO in the Schrodinger equation of the form suggested by 
BenDaniel and Duke (/3 = — 1, 7 = 0) 0, we conclude that both equations lead to the same 
solutions in the energy range concerned. 

Let us now introduce the numerical method that allows us to obtain the energy levels 
for a Dirac equation with space-dependent effective mass. 

We will consider a particle with mass m{z) that is submitted to an arbitrary one dimen- 
sional potential well V{z). The time-independent Dirac equation is written as (in units with 

n = c = i)m 



{ap + I3m{z))^ = {E -V{z))^, (2) 

is the m( 

matrices given by 



where, p = —i-r is the momentum operator, E is the electron energy, a and (5 are 4x4 

dz 



l^a^ j' ^=1,0 -/J' 

/ is the 2x2 identity matrix and cr^ is a 2 x 2 Pauli matrix defined as 



-^ = 1;:) (4) 




Consider now an arbitrary well as sketched in fig.l. We split up the interval [a, 6] into 
infinitesimal intervals of length Az = {b — a)/N. For the i-th interval, we approximate 
the potential and the mass by 

V{z) = V{zi) = Vi and m{z) = m{zi) = mi, for Zi < z < Zi^i. (5) 

The wave function of the electron in Dirac's equation with no spin flip in the i—th interval 
is 
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where pi = \J {E — Vj) — mf. The bound state conditions are given by \E — Vq\ < mo and 
\E — Vat+iI < "n^AT+i- By imposing the continuity of the wave function at each z = Zi, we 
have a matrix M{E) which relates the coefficients in the region where z < a with the region 
where z > b. 




The finiteness of the wave function requires that 

M{E)22 = 



(7) 
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So, the solution of (||) give us the energy levels. 

Note that this numerical method is specially convenient for treating wells and barriers 
with arbitrary profiles and it is nothing else than an extension of the transfer-matrix method 
for relativistic theories. To the best of our knowledge, this is the first numerical analysis 
for evaluating energy levels in an relativistic equation with mass position-dependent and 
arbitrary potential. However, our main point here is to demonstrate that the Dirac equation 
can be used to obtain unambiguously results in situations where the Schrodinger equation 
depends on ordering problems. 

As an illustration, we applied the method described above for an electron in a one- 
dimensional GaAs/Alo.sGao.rAs heterostructure. For the sake of comparison with previous 
results we take a square well, as sketched in fig. 2. The electron effective mass is 0.67mo and 
0.86mo for GaAs and Alo.3Gao.7As respectively ( niQ is the free electron mass) (Fig. 2). In 
the Schrodinger equation we use the von Ross operator considering several values of the a 
parameter. Note that for abrupt heterojunctions only Hamiltonians with = 7 are viable, 
due to continuity conditions across the heteroj unction 0. As we can see in fig. 3 there is 
an extraordinary coincidence between the results from the Dirac and from the Schrodinger 
equations for a = 0. As a matter of fact, this is expected since the maximum value of the 
energy involved is 3eV. Further, this result strongly supports the prescription of BenDaniel 
and Duke |^ for the Schrodinger context. 

As a second illustration consider the conduction-band structure of a GaAs/Alo.aGao.yAs 
system with nonabrupt interface and assume that the effective mass changes linearly at the 
transition regions while the potential well varies quadratically in that regions (denoted by 
a) as it is shown in fig. 4. There the potential is given by 



Viz) = C ea(^)+e2X(^)' 



(9) 



where C = 0.6 is the conduction band offset and ei = 0.3, €2 = 0.7 are constants associated 
with the compositional dependence of the energy-gap difference between GaAs and AlGaAs 
(experimental parameters and details concerned can be seen in refs. f^J^l)- 
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Again, we use our relativistic method and the Schrodinger equation with the BenDaniel 
and Duke prescription (a = 7 = 0). Once more, as we can see in fig. 5, a complete 
coincidence between the relativistic and non-relativistic results is obtained. 

In conclusion, we have shown that a relativistic method can be successfully used to 
overcome the ordering problem of the kinetic energy operator in non-relativistic models. 
Since the range of energy involved is extremely low (comparing to electron rest mass), the 
numerical results are perfectly coincident in both cases. 

It is worthwhile to mention that, notwithstanding the Wannier-Slater theorem com- 
mented in the introduction, we claim attention for the coincidence shown above. Therefore, 
we believe that, after exhaustive and appropriate considerations about effective-mass approx- 
imation, band structures and the periodic potential, the relativistic approach constructed 
here it could be used to calculate physical parameters in the theory of abrupt and nonabrupt 
semiconductor heterostructures. 

Moreover, the relativistic treatment developed here can be applied to all physical systems 
described by a Sturm-Liouville eigenvalue equation (within an appropriate range of energy), 
namely 



d 

dz 



' 1 df 
p{z) dz 



+ q{z)f\z) = \w{z)f{z) . (10) 



Here / is an eigenf unction, A is an eigenvalue and p, g, and w describe particular properties 
of the system. For example, this equation can describe the motion of electrons or phonons 
along the growth axis of a [100] zinc-blende heterostructure |jlO|. Thereby, the equation (p^ 



can always be replaced by an correspondent Dirac equation for avoid ordering problems. 
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FIGURE CAPTIONS 



Figure 1: Generic profile of a one- dimensional quantum potential well. We split up the 
interval [a, b] into N infinitesimal intervals of length Az — {b — a)/N 

Figure 2: Potential well and effective mass, as a function of the position, for an abrupt 
heterostructure. 

Figure 3: Eigenenergies in the conduction band of the GaAs/Alo.sGao.rAs, with a con- 
duction band-offset of 0.6 vs the well width (abrupt hetero junction). The chosen value for 
a is 20% of the well width. Solid line shows the calculations performed using the Dirac 
equation. Dotted line shows the calculations for Schrodinger equation using the BenDaniel- 
Duke prescription (a = 0) . The -|- -|- curve denotes calculations using the Zhu and Kroemer 
prescription (a = —0.5). 

Figure 4: Potential well and effective mass, as a function of the position, for a nonabrupt 
heterostructure ( a denotes the transition region). 

Figure 5: Eigenenergies in the conduction band of the GaAs/Alo.sGao.yAs, with a conduc- 
tion band-offset of 0.6 vs the well width (nonabrupt hetero junction). The chosen value for a 
is 20% of the well width. Solid line shows the calculations performed using the Dirac equation 
and dotted line shows the calculations for Schrodinger equation using the BenDaniel-Duke 
prescription (a = 7 = 0). The -|- -|- curve denotes calculations using the Zhu and Kroemer 
prescription (q; = 7 = — 0.5). 
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